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CN| ABSTRACT 

$_ i In this paper, the detectability of habitable exomoons orbiting around giant planets 

Qh m M-dwarf systems using Transit Timing Variations (TTVs) and Transit Timing Du- 

«^ rations (TDVs) with Kepler-class photometry is investigated. Light curves of systems 

with various configurations were simulated around M-dwarf hosts of mass 0.5 M and 
radius 0.55 Rq. Jupiter-like giant planets which offer the best potential for hosting 
habitable exomoons were considered with rocky super-Earth-mass moons. The de- 

i— | tect ability is measured by using the phase-correlation between TTV and TDV signals. 

Qn Since the TDV signal is typically weaker than the TTV signal, confirmation of an 

nQ exomoon detection will depend on being able to detect a TDV signal. We find that 

exomoons around planets orbiting within the habitable zone of an M-dwarf host star 
can produce both detectable TTV and TDV signatures with Kepler-class photometry. 
While aliasing between the planet period and moon period may hinder exomoon detec- 
tion, we also find some strong correlation signatures in our simulation (eg. correlation: 
>0.7) which would provide convincing exomoon signatures. With the addition of red 



Oh 



$H 



c/} noise stellar variability, correlations generally weaken. However simulated examples 

, ^, with planet masses less than around 25 M@, moons of mass 8-10 M e and specific 

values of planet and moon periods still yield detectable correlation in 25-50% of cases. 
^H Our simulation indicates that Kepler provides one of the best available opportunities 

^ for exomoon detection. 
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CO 1 INTRODUCTION timing (Lewis et al. 2008), Rossiter-McLaughlin effect (Si- 

^H mon et al. 2010) and scatter-peak (Simon et al. 2012). 

• • Over the last decade, the search lor and study or exopianets A . . r ', ., ,, , . , . , ., , . 

>_ ' . J An extension or the transit method involving transit tim- 

is one ol the most dynamic research fields ol modern astron- . . rr ,, , , , . , , , , , , , ., , , 

. _ _ _ . ;:_.,_ . ™ i ! , ing appears to oner the best potentral to detect habrtable 

k> omy. As ol March 2013, more than 850 planets have been . ,, - , ,,,. . nni1l N ™ 

J\] J ill! i i i exomoons in the near luture (Kipping 201 lb). Ihe pres- 

r , ^ confirmed Irom ground-based observatron and more than - . , . , . . ,, 

_ „__,_ . .. . . . .. . . . , _ . ence ol an exomoon induces two mam variations on to the 

rn 2,300 planet candidates have been discovered by the Kepler , , , , m -*. m- • tt • j.- /mmT 7 \ /o 

^ . . /rk „ . , A10 , . . i r i i host planet; Transit Timing Variation (TTV) (Sartoretti 

mission Batalha et al. 2013). As the number ol detected a , . , in nn\ j r^ •. t^ .- Tr • \. /m^A 

_ v ' ii-i & Schneider 1999) and Iransit Duratron Varratron (1DV) 

exopianets continues to grow, the potential lor detecting /T ,. . ™™ i \ ^^t j r^™r i v ^ w 

,.. . . . . ii n • • (Kipping 2009a, b). TTV and TDV signals are predicted to 

satellites orbiting them has become ol increasing interest , /0 , £ , r , , , ., 

° . _ . . . -~~~ TT . TT be 7T/2 out ol phase lor an edge-on circular orbit system, 

in recent years bartoretti & Schneider 1999; Han & Han ' . . T , , , . ,, 

^ r v , ,., , x ' creating a unique exomoon signature. In order to claim the 

2002; Kipping 2009a,b; Kipping et al. 2009). The presence r ., 6 , w . . , , . , ,, 

_ . i i i -t r i • presence ol an exomoon with the timing technique, both 

ol exomoons may improve the probability ol the existence mrn^r i mnw i , t i . , i 

-.._ I./ i ii i i i TTV and TDV signal must be detected, 

ol hie on their host planet and the moons themselves also 

have potential to host life (Laskar et al. 1993). Moreover, the Ideal host stars for habitable exomoon detection are M- 

detection of exomoons would improve our understanding of type stars, due to the large amplitude of transits generated 

planetary formation and evolution (Williams et al. 1997). (Charbonneau et al. 2009; Bean et al. 2010; Vogt et al. 2010; 

Several methods have been developed to detect exo- Mann et aL 2012 ) and the sma11 distance of their habitable 

moons; including the transit (Simon et al. 2007), microlens- zone ; which increases the transit probability and the number 

ing (Han & Han 2002; Liebig & Wambsganss 2010), plusar- of transit events per observation time (Gaidos et al. 2007; 

Kaltenegger 2010). Several studies (eg. Barnes & O'Brien 

(2002); Domingos et al. (2006)) explore the stability of orbits 

* E-mail: supachai.awiphan@students.manchester.ac.uk around gas giants. In 2002, Barnes & O'Brien proposed that 

1 See http://exoplanet.eu/ Earth-mass moons of Jupiter-mass planets around stars of 
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mass greater than 0.4 Mq are dynamically stable for billions 
of years, long enough to sustain life. 

At present, the best instrument up to this challenge is 
Kepler which launched in 2009. Kepler is designed to detect 
small transiting exoplanets with its highly sensitive photo- 
metric camera. It is monitoring ^150,000 stars in Cygnus, 
including more than 3,000 M-dwarf stars (Batalha et al. 
2010). Many studies in exomoon detection with Kepler have 
been done (eg. Szabo et al. (2006); Kipping et al. (2012); Si- 
mon et al. (2012)), including a study of Kipping et al. (2009) 
using TTV and TDV techniques that found that exomoons 
around gas giants planets in the habitable zone of M-dwarf 
stars may be detectable with the Kepler Mission. The Hunt 
for Exomoons with Kepler (HEK) project aims to search 
Kepler data for evidence of exomoons (Kipping et al. 2012). 

In this paper, we assess the detect ability of habitable 
exomoons using the phase-correlation between TTV and 
TDV signals from simulation of Kepler-class photometry. 
The structure of the paper is as follows. In Section 2, a 
background of habitable exomoons with their timing effects 
is provided. In Section 3, the target selection is described. 
The generating of light curve of an exoplanet with an exo- 
moon with Kepler is described in Section 4. The last part 
of this Section is dedicated to analysing the detectability of 
habitable exomoons using a detailed numerical simulation. 
In Section 5, simulation results are presented and analysed. 
Finally, the methods and results are summarised in Section 
6. 



2 HABITABLE EXOMOON DETECTION 

2.1 Habitable exomoons 

An exomoon is a natural satellite of an exoplanet. None 
have been discovered to date. However, if our Solar System 
is typical, then exomoons must be common. 

In 1997, Chyba and Williams et al. proposed the pos- 
sibility of habitable exomoons. They proposed that moons 
orbiting around giant gas planets' Hill sphere could host life. 
Although moons orbiting giant planets at 1 AU from a solar 
analogue would become tidally locked within a few billion 
years after they form, their orbital period on timescales of a 
few days to a few months could cause temperature fluctua- 
tions on them (Chyba 1997; Williams et al. 1997). Kipping 
et al. (2009) suggested that the habitable zone of exomoons 
can be defined simply as the distance where planets receive 
the same energy as the Earth, 



magnetic activity. Therefore, if planets or moons around the 
dwarfs can hold on to their atmosphere in the early age and 
are not tidally locked, they could be habitable (Joshi et al. 
1997). 

Habitable moons need to be large enough to retain wa- 
ter and an atmosphere. Although moons formed from the 
proto-planetary disk are unlikely to be greater than 0.02% 
the mass of the host planet (Canup & Ward 2006), moons 
formed from captures (Triton; Agnor & Hamilton (2006)) 
or impacts (the Moon; Taylor (1992)) have no limit on their 
mass. The results of Weidner & Home (2010) show that 
there are some planets which could retain massive moons 
(Earth-size moon), including detected exoplanets such as 
CoRoT-3b and CoRoT-9b. 

Another major factor that should be taken into account 
is orbital stability. Because of three-body instability, moons 
can be lost from their host planet if the distance between the 
planet and moon is too large. To be retained by the planet, 
the moon must have an orbit that lies within the Hill sphere, 
Rh (Barnes & O'Brien 2002). 



Rh — CLp( 



M p 
3M* 



1/3 



(2) 



clhz = 



AU. 



(1) 



where M* is mass of star, M p is mass of planet and a p is star- 
planet semi-major axis. In this project, the satellites within 
the Hill sphere of the habitable zone planets are considered, 
due to their habitability. 

In the remainder of this paper, we use the terms "moon" 
and " exomoon" to refer to the less massive companion of a 
planet-satellite system even where the companion may itself 
be of planetary mass. 



2.2 Transit Timing Variations (TTV) 

In order to detect exomoons, several methods have been de- 
veloped. The TTV and TDV techniques are focused upon 
in this work since they could detect small exomoons (Sar- 
toretti & Schneider 1999; Kipping 2009a,b). The concept of 
the TTV technique is that the presence of a third body such 
as exomoon in the system causes a change in planetary or- 
bit. The time between transits varies because the transiting 
planet and the moon exchange energy and angular momen- 
tum. This gravitational interaction perturbs the orbit of the 
transiting planet and causes a short-period oscillation of the 
semi-major axes and eccentricities. The signal depends on 
the mass, separation and orbital parameters of the planet 
and the moon (See Kipping (2011b)). For edge-on circular 
orbits, the TTV signal and RMS amplitude of the signal can 
be written as, 



where L* is the host star luminosity. However, stellar evolu- 
tion can alter the boundaries of the habitable zone. Due to 
stellar evolution, the host star becomes brighter and hotter 
which shifts the habitable zone outward (Selsis et al. 2007). 
For high-mass stars, their habitable zones are much further 
and broader which could be the best candidates for find- 
ing habitable planets or moons. However, emission in the 
far ultraviolet (FUV) band which is potentially damaging 
to life is observed in high-mass stars. While low-mass stars 
have the longest lifetimes, in their early age they have strong 
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where M m is the moon's mass and a m is the semi-major axis 
of the moon around the planet-moon barycentre, P p is the 
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orbital period of the planet and f m is the true anomaly of 
the moon. 

Unfortunately, TTV can also be induced by a multitude 
of phenomena, including general relativistic precession rate 
of periastra (Jordan & Bakos 2008), stellar proper motion 
(Rafikov 2009) and parallax (Scharf 2007) effects. Therefore, 
a TTV signal by itself cannot confirm the presence of an 
exomoon. 



2.3 Transit Duration Variation (TDV) 

TDV is the periodic change in the transit duration over 
many measurements caused by the apparent velocity of the 
planet which increases and decreases due to the planet-moon 
interaction. Kipping showed that exomoons should induce 
not only the TTV effect but also the TDV effect on their 
host planets. 

For the systems with non-coplanar orbits, the TDV ef- 
fect can be separated into two main constituents, a velocity 
(V) component and a transit impact parameter (TIP) com- 
ponent (Kipping 2009b). The V-component is caused by the 
variation in velocity of the planet due to the moon's grav- 
ity. The TIP-component is affected by the planet moving 
between high and low host-star impact parameters. Kipping 
(2011b) formulated that the total TDV signal as a linear 
combination between TDV-V signal and TDV- TIP signals 
(See Kipping (2011b)). 

In a system with 90 degree orbital inclination and a 
circular orbit, the TDV TIP-component will be zero. There- 
fore, the TDV signal of edge-on circular orbit is, 



TDV 



Q>m-Lvl m Tp 



CLpMpPrr 



sin(/ rr 



(5) 



well as the transit frequency per observation time (Kalteneg- 
ger 2010). Sasaki et al. (2012) also suggested that the semi- 
major axis of the host planet for the most detectable ex- 
omoons around an M-dwarf star is 0.2-0.4 AU. Therefore, 
the most detectable exomoons in M-dwarf systems can or- 
bit within the habitable zone. However, the host with mass 
less than 0.2 Me cannot host a habitable exomoon (Heller 
2012). 

Kepler monitored the Cygnus region along the Orion 
arm centred where there are about 0.5 million stars brighter 
than 16 th magnitude (Kepler passband) within its FOV. 
However, only 10 5 stars with magnitude less than 16 are ex- 
pected to be exoplanet hosts. In 2010 the Kepler mission an- 
nounced 150,000 highest priority target stars, but only 2% of 
these target stars have effective temperature less than 3500 
K (Batalha et al. 2010), whereas >70% of all stars within 
20 pc are M-dwarfs (Henry et al. 1994; Chabrier 2003; Reid 
et al. 2004). However, in 2011, the team released additional 
exoplanet data, including 997 planet-candidate host stars in 
which 74 (>5%) have effective temperature less than 4400 
K in the Kepler Input Catalog (Batalha et al. 2010; Borucki 
et al. 2011; Brown et al. 2011). 

For our simulation of TTV and TDV signals, we as- 
sume that the host is an M-dwarf star with mass 0.5 Mq 
and radius 0.55 Rq. Their effective temperature, microtur- 
bulent velocity and logg are set to be 3500 K, 1 km.s -1 and 
4.5, respectively, as applicable to solar-metallicity M-dwarf 
(Bean et al. 2006; Onehag et al. 2012). In order to calculate 
the limb darkening coefficient, solar-metallicity is assumed 
and a quadratic limb- darkening model is used. The values 
of limb- darkening coefficients for the transmission curves of 
Kepler are obtained from Claret & Bloemen (2011). For our 
M-dwarf targets, the value of the coefficients 71 and 72 are 
0.4042 and 0.3268, respectively 2 . 



and the RMS amplitude of the TDV signal is, 



_ ClrYilVlfYiTp 

CLpMpPmVZTT 



(6) 



where P m is the orbital period of the moon and f is transit 
duration. The TDV technique cannot detect habitable ex- 
omoons alone because the TDV signals are relatively weak 
compared with the TTV signals (Porter & Grundy 2011) 
and can also be induced by parallax effects (Scharf 2007). 
However, combining TDV and TTV signals can confirm the 
presence of an exomoon, because the signals have a ^ -phase 
difference that provides a unique exomoon signature. The 
orbital separation and mass of exomoons can also be ob- 
tained. 



3 MODELLING HABITABLE EXOMOONS 



3.2 Properties of the host planet 

Jupiter-like giant planets offer the best potential for detect- 
ing habitable exomoons (Kipping et al. 2009). In order to 
investigate habitable exoplanets and exomoons, the planet- 
star separation is set to be inside the habitable zone, starting 
at a separation of 0.10 AU and increasing in logarithmically 
to 0.66 AU. This range includes the semimajor axis of M- 
dwarf planets which Sasaki et al. (2012) argue on stability 
ground may be among first detectable exomoon systems. 

We simulate giant planets with masses ranging loga- 
rithmically from 15 to 150 M©. Fortney et al. (2007b) found 
that the radius of giant planets depends on their overall 
mass, core mass and separation. For giant planets of age 
4.5 Gyr, their radius falls between 1.0 and 1.2 Rj (Jupiter 
radius). Therefore, we adopt a planet radius of 1.2 Rj. 



3.1 Properties of the host star 

In this analysis, M-dwarf stars are selected to be the exo- 
planet host stars. Very cool (late K and early M type) dwarf 
stars have become popular targets of planet searches, be- 
cause the amplitudes of the transits generated by planets 
in M-dwarfs are larger than those generated by hotter stars 
(Charbonneau et al. 2009; Bean et al. 2010; Vogt et al. 2010; 
Mann et al. 2012) and the small distance of their habitable 
zone increases the transit probability of habitable planets as 



3.3 Properties of the exomoon 

No exomoon has yet been discovered, therefore the proper- 
ties of Earth-like planets are used for the habitable exomoon 
in this work. Rocky planets with logarithmic mass between 



See http : //cdsarc . u-strasbg . f r/viz-bin/qcat? J/A+A/529/A75 
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1 and 10 Mq are chosen. The radius of the moon is calcu- 
lated from Fortney's model, using rock mass fraction equal 
to 0.66 (Earth-like planet): 



Table 1. Keller photometry properties. 



Rm = 1.00 + 0.65 log M m + 0.14(log M m ) 



(7) 



where M m and R m are the moon's mass and moon's radius 
in Me and R®, respectively (Fortney et al. 2007a,b). Only a 
moon within the planet's Hill sphere with an orbital period 
between 1.00 to 3.16 days is considered. Again, for simplicity, 
circular orbits are assumed. 



4 GENERATING AND ANALYSING LIGHT 
CURVES 

4.1 Kepler transiting light curve generation 

In order to generate transit light curves of a planet with a 
moon, the algorithms of Kipping (2011a) are used. The Ke- 
pler mission is designed to monitor ^150,000 stars brighter 
than 16 t/l magnitude (in the Kepler passband) with 20 parts 
per million photometric precision at 12 t/l magnitude in 6.5 
hours (Batalha et al. 2010; Caldwell et al. 2010b). In order 
to meet this requirement, the estimated photon collection 
rate of Kepler is (Borucki et al. 2005; Yee & Gaudi 2008; 
Kipping et al. 2009), 



- ph 



6.3 x 10 s fi 10 



-0.4(m-12) 



photons/hours, (8) 



where m is the apparent magnitude. However, Kepler pho- 
tometry is also affected by shot noise, background flux and 
instrumental noise. Table 1 summarise the properties we as- 
sume for Kepler photometry, including noise contributions 
we now discuss. 



4-1-1 Shot noise 

Shot noise or Poisson noise comes from the discrete nature 
of photons. At 12 t/l magnitude, the largest noise component 
is the Poisson noise of the target (Caldwell et al. 2010b) 
which we simulate. 



4-1-2 Background flux 

The background flux for Kepler comes from zodiacal light 
from the Solar System and diffuse starlight from background 
stars. In pre-launch prediction, the background flux is esti- 
mated at around 334 e~sec _1 or 22 magnitudes per square 
arcsecond (Caldwell et al. 2010b). However, in real observa- 
tions, the background flux varies across detectors and with 
orientation of the telescope. We adopt the pre-launch back- 
ground flux estimate to generate the light curves in this 
work. 



4-1-3 Instrumental noise 

There are two main components of instrumental noises for 
Kepler: read noise and dark current. From in-flight measure- 
ment, the read noise median value is 95 e~read _1 and dark 
current is 0.25 e~pixel _1 s _1 which is quite low compared 
to the photons collected from the targets (Caldwell et al. 



Parameter 


Value 


Exposure time (s) 


6.02 


Plate scale (arcseconds/pixel) 


3.98t 


Background flux (e _ s _1 ) 


334t 


Read Noise (e _ read _1 ) 


95* 


Dark Current (e — pixel - 1 s —1 ) 


0.25* 


* : Van Cleve & Caldwell (2009) 




* : Caldwell et al. (2010a) 




t : Caldwell et al. (2010b) 





2010a). We include it in our simulation, though its affect on 
our results is negligible. 



4.2 Measuring TTV-TDV signals 

In order to find the transit time of minimum and transit du- 
ration, the ingress and egress of our simulated light curves 
are fitted. The light curves are divided into phase bins using 
the input period which is assumed to be precisely deter- 
mine from observational data. A running straight-line fit is 
made to three consecutive points of phased data. The fits 
with minimum and maximum slopes are chosen to define 
the ingress and egress of the transit, respectively. The inter- 
section between the light curve median and the ingress and 
egress slopes are used to define the ingress (ti ng ) and egress 
(tegr) times, respectively (Figure 1). The time of minimum 
light (to) and the transit duration (f mea n) are defined as, 



to = (ting + t e ^ r )/2, and Trr. 



respectively. 



Using the mid-transit time, a new ephemeris as a func- 
tion of epoch is derived. The new ephemeris is determined 
by fitting a linear function to the mid-transit points. 



To(n)=T o (0)+nP 



(9) 



where n is epoch and To is time of minimum light as a func- 
tion of epoch. The residuals of the times of minimum and 
transit duration taken as the TTV and TDV signal of the 
system. 



4.3 TTV-TDV correlation testing 

Holman & Murray (2005) showed that other planets could 
induce TTV signals on a transiting planet. Therefore, the 
TTV signal alone cannot distinguish between the effect of 
other planets and the effect of exomoon. In order to confirm 
an exomoon TDV signals must also be detected. 

From Section 2.3, TTV and TDV signals are sinusoidal 
functions and the TTV signal is 90 degrees out of phase with 
the TDV signal in coplanar systems. In the case of a circular 
planetary orbit and a co- aligned moon orbit, the TDV-TIP 
component exists. From Equation 3 and Equation 5, the 
TTV signal, TDV signal and the relation between TTV and 
TDV are, 



TTV = 



OlmlVlmip 

27ra p M p j 



COS (frr 



(10) 
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Figure 1. Light curves of a 15 ¥9 habitable-zone planet with a 10 ¥9 moon for an M2 host star with planet period 89.35 days 
and moon period 2.24 days with differing moon phase. Error bars are shown at 1,000 times their true size. The fits show ingress 
time and egress time of transit, and average flux (median of flux data). 



TDV=f 



and 



TDV = 
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(11) 



dv-Lrr, 






+ 



1-6? R, 



l-bl R* 



TTV 



(12) 



where b p is impact parameter of the planet from its host and 
R* is host star radius. However, the ratio between TDV-V 
and TDV- TIP is very large. For the systems in our simula- 
tion, the minimum ratio is 1,100. Therefore, the TDV- TIP 
is negligible and the relation between TTV and TDV signal 
can be written as, 



TDV 2 



/gTTfV 



- I TTV 2 + f 2 

\ ± m J 



Q>m-Lvl m ±p 
dp MpP m 



(13) 



Therefore, in theory, the plot between the square of the 
TTV signal and square of the TDV signal should show a per- 
fect linear relationship with negative slope. However, there 
are other effects, such as star spots, instrument noise and 
sparsity of observation which could produce false positive 
TTV and TDV signatures. 

In this simulation, the instrument noise and observing 
frequency both affect the TTV and TDV signals. Thus, the 
plot between TTV 2 and TDV 2 may not show a clear linear 
relationship. In order to check this relationship, the Pear- 
son product-moment correlation coefficient was calculated 
to test the correlation between TTV 2 and TDV 2 . The co- 
efficient is 



X = 



J27=i(TTV? - TTV 2 )(TDV 2 - TDV 2 ) 



J2?=i(TTV 2 - TTV 2 ) 2 ^/J2?=i(TDV 2 - TDV 2 ) 2 

(14) 



A negative coefficient is produced by an inverse rela- 
tionship between the two variables and a positive coefficient 
means there is a positive linear relationship. The positive 
slope of TTV 2 and TDV 2 plot means the TTV and TDV 
signal are not consistent with sinusoidal functions with a 90 
degree phase difference. The TTV signal, TDV signal and 
TTV 2 versus TDV 2 of three sample systems are shown in 
Figure 2. 



5 DETECTABILITY OF HABITABLE 
EXOMOONS 

5.1 Detectability of habitable exomoons 

The light curves are generated with Kepler photometric 
noise. 146,410 light curves are simulated with 11 indepen- 
dent values of each of four variable input parameters: planet 
mass; planet separation; moon mass; and moon period, and 
10 random initial phases. The host stars are assumed to be 
M-dwarf stars of 12.5 magnitude in the Kepler passband. 
The cadence of this simulation is 50 data points per day 
(every 28.8 mins) which corresponds closely to Kepler's long 
cadence mode (every 29.4 mins) (Gilliland et al. 2010). In 
order to simulate the current Kepler data, a 3-year simu- 
lation of a transiting giant extrasolar planet with a rocky 
extrasolar moon was run to find out the detectability of 
an exomoon in the M-dwarf habitable zone. The details of 
physical parameters of the systems are listed in Table 2. 

The 4D-simulation is projected on to two-parameter 
planes in order to examine the relation between two vari- 
ables. Since we are only interesting in negative correlations, 
we define the projected correlation as: 



Xproj 



1 

N 



5Z Xi ' 

i( X >0) 



(15) 



where N is the total number of 2-D simulations that are 
projected and i(x > 0) refers only to those simulations with 
negative correlation. The projected plots therefore represent 
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Figure 2. TTV signal (left), TDV signal (middle) and TTV 2 versus TDV 2 (right) for a 10.0 M e exomoon of a planet orbiting 
around a 0.5 solar-mass M-dwarf star. The with planet period is 89.35 days and moon period is 2.24 days. Panels are shown in 
three exoplanet with log(M p /15M e ) equal to 0.0 (a), 0.4 (b), and 0.7 (c). 



averages over logarithmic parameter priors for negative cor- 
relation signals. In the left hand panels of Figure 3(a), the 
plot between planet mass and moon mass shows that a high- 
mass moon hosted by a low-mass planet is the most de- 
tectable of the systems considered. This result agrees with 
the moon period versus planet mass and moon mass plots 
(Figure 3(e) and Figure 3(f)). However, in these two plots, 
the changes in moon period do not affect the correlation. 
In Figure 3(b), the projection plot between separation of 
planet and period of moon also does not show any signifi- 
cant trends. 

Figure 3(c) shows the detectability coefficient between 
mass and separation of the planet. Planets with high sepa- 
ration have higher detectability than close-in planets of the 
same mass. This result correlates with the result of moon 
mass versus planet separation (Figure 3(d)) which shows 
that, in systems of equal satellite mass, the outer planet 
hosts have larger correlation coefficients. These features may 
be produced by only a few transit events in high planet sep- 
aration systems, because, at 0.6 AU separations, only three 



transit events are detected in the simulation. Therefore, we 
now check the reliability of the correlation. 

The analysis of correlations is meaningful when the cor- 
relations are not dominated by noise. The variance of cor- 
relation is plotted in the right-hand panels of Figure 3 in 
order to check the reliability of testing. From Figure 3(a) to 
Figure 3(d), the variance plots show that the systems with a 
small number of transit events (long planet period systems) 
have higher variance. However, the value of the variance is 
still low compared to the correlation coefficient. 

While the magnitude of Xproj in Equation 15 is reduced 
by positive correlations that are included in A", we have 
checked that the basic features in the plots of Figure 3 trace 
those obtained by ignoring positive correlations, albeit at a 
weaker level. Finally, the assumption that moons of outer 
planets should be easier to detect than moons of inner plan- 
ets is confirmed by our simulations. 

The theoretical lines of RMS amplitude of the TTV 
and TDV signals are shown to investigate the features of 
the contours. For Figure 3(a) the high amplitude of TTV 
and TDV signals produces a high coefficient of detection 
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Table 2. Input parameters assumed for our exomoon simulations. 



Star parameters 


Mass (M©) 


0.5 


Radius (Rq) 


0.55 


Apparent magnitude (K p ) 


12.5 


Quadratic limb-darkening coefficient 1 


0.4042 


Quadratic limb-darkening coefficient 2 


0.3268 


Planet parameters 


Mass (M e ) 


15.0-150.0 


Radius (Rj) 


1.2 


Separation (AU) 


0.10-0.66 


Eccentricity 


0.0 


Inclination (degrees) 


90.0 


Moon parameters 


Mass (M e ) 


1.0-10.0 


Radius (-R©) 


Equation 7 


Period (days) 


1.00-3.16 


Eccentricity 


0.0 


Inclination (degrees) 


90.0 



with the same slope. Moreover, the features in Figure 3(c) 
and Figure 3(d) are also well- correlated with TDV RMS 
amplitude signals which can be explained by the relative 
weakness of TDV signals compared with TTV signals. In 
conclusion, the detectability of exomoons is dominated by 
the amplitude of TDV signals. 



5.2 Analysing the correlation structure 

The structures of the correlation plots are explained by the 
magnitude of the TDV signal. However, in Figure 3(c) and 
Figure 3(d), gaps are evident at planet semi-major axes of 
0.4 and 0.5 AU. The variance plots show that there is no 
difference in variance across this region and that therefore 
the features in Figure 3(c) and Figure 3(d) are real. To 
investigate the structures, the 4-D plots are sliced into 2- 
D plots. The correlation plots of planet separation versus 
planet mass, and of planet separation versus moon mass are 
shown in Figure 4 and Figure 5. Systems with high moon 
mass and low planet mass have a high value of correlation 
and the features in these contours correspond to the pro- 
jected contours of Figure 3, including the gap structures. 

In Figure 4 and Figure 5, the maps with similar moon 
period show gap features at the same planet separation, 
but they shift with a different moon period. The ratios be- 
tween moon period and planet which produce the gap are 
near-integer values and correspond to the cold spots in Fig- 
ure 3(b). Therefore, the gap structures can be explained by a 
resonance between the planet and moon period which pro- 
duces constant detected TTV and TDV signals. However, 
they also depend on the number of detected transit events. 
In short period systems which have a larger number of tran- 
sits, the gap structures are more difficult to produce due to 
the larger range of detected planetary phases. 



5.3 M-dwarf variability 

In the previous section, the light curves without intrin- 
sic stellar variability were simulated. In reality, stars can 
vary in in brightness due to pulsation, rotation and activ- 
ity (red noise). Around 40-70% of M-dwarfs have variability 
with photometric dispersion (a m ) ^3-5 mmag, depending 
on their brightness. For stars with 12.5 magnitude in the 
Kepler passband, the variability fraction is nearly 1 (Ciardi 
et al. 2011) and their noise tends to have long variation pe- 
riods (^5 days) (McQuillan et al. 2012). Therefore, in our 
simulation, the red noise with a 12 day period was added 
into the light curves in order to investigate the effect of stel- 
lar noise to the detectability. Their amplitude based on that 
found for M-dwarfs with Kepler magnitude between 12-14 
(Ciardi et al. 2011). In order to simulate short term vari- 
ability, five minor variations with linearly random periods 
between and 12 days and amplitudes less than half of the 
amplitude of main variation were also added. An example 
of our simulated red noise is shown in Figure 6. 

Our test of the effects of red noise were based on follow- 
ing parameters. We simulated planets with log(Mp/15Me) 
equal to 0.0, 0.1, 0.2, 0.3, 0.4 and 0.5 and tested moons with 
log(M m /M©) equal to 0.8, 0.9 and 1.0 M®. The planet and 
moon are assumed to have a period of around 89 and 2.2 
days, respectively, that correspond to the peak in Figure 4 
and Figure 5. For each parameter combination, we used sim- 
ulations with high TTV-TDV correlation (>0.7) to exam- 
ine the effect of red noise on high-confidence detections. We 
added 500 different variations for each dispersion to each 
light curve. The result in Figure 7 shows that the presence 
of intrinsic stellar variability of M-dwarfs might effect the 
exomoon detectability. The stellar variability reduces the 
exomoon detection correlation by 0.1. However, for our sim- 
ulated systems with planet masses less than around 25 M e 
with moon masses 8-10 M®, typically 25-50% of them still 
have correlations high enough to confirmed exomoon detec- 
tion. 



5.4 Improving the detectability with long-term 
observation 

Section 5.1 and 5.2 are based on 1,000-day simulations of 
50,000 data points (every 28.8 minutes), which corresponds 
to the integration time of the Kepler long-cadence strategy. 
In theory a longer time baseline should produce a better 
detectability coefficient, because of more transiting events. 
In 2012, Kepler has been approved for extension through 
to 2016 . Therefore, 5-years photometric data from Kepler 
targets should be obtained. We have performed a 1,600-day 
simulations to produce a sample of long-term observations. 
The moon period is set to be 2.51 days, because it is not 
in resonance with any planetary period considered in this 
work. 

The contours of two simulations show the same struc- 
tures in Figure 8, but the long-term simulation has a better 
detectability. From this result, the number of detected exo- 
moons should increase with length of observation. The ex- 
tension to the Kepler mission and possible additional space- 
based telescope missions (eg. PLATO (Rauer 2012)) should 
allow the habitable exomoons to be detected or their abun- 
dance constrained. 
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(a) Planet mass versus moon mass 
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(b) Planet semimajor axis versus moon period 
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(c) Planet mass versus planet semimajor axis 



Figure 3. Correlation (left) and variance of correlation (right) between (a) planet mass and moon mass, (b) planet semimajor 
axis and moon period, (c) planet mass and planet semimajor axis, (d) moon mass and planet semimajor axis (next page), (e) 
planet mass versus moon period, and (f) moon mass and moon period of the light curves. The contour is averaged over other 
two variable. The RMS amplitude of the TTV signal (dashed) and RMS amplitude of the TDV signal (dot-dashed) in units of 
seconds are presented. The cold spots in (b) indicate the data with the planet period in resonance with the moon period (See 
Section 5.2). 
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(d) Moon mass versus planet semimajor axis 
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(e) Planet mass versus moon period 
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(f) Moon mass versus moon period 
Figure 3. Continued 
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Figure 6. Our simulated stellar red noise with a main noise com- 
ponent of 12 days period (See Section 5.3). 



6 CONCLUSION 

In this work, the light curves of a transiting exoplanet with 
an exomoon were implemented for the purpose of determin- 
ing detectability of exomoons. The Keller photometric noise 
was modelled to the light curve in order to simulate the data 
from Kepler. Measuring the detectability was done by phase- 
correlation between TTV and TDV signals. TTV and TDV 
always exhibit a 90- degree phase shift, therefore, the TTV 2 
signal is linear with the TDV 2 signal. The Pearson product- 
moment correlation coefficient was used to determine the 
detectability of signals. 



3-year Kepler light curves of 146,410 systems with 
various configurations were simulated. For each extrasolar 
planet system, the giant planets and their rocky satellites 
were placed in the habitable zone of 0.5 Mq dwarf stars. 
Their masses and periods were selected logarithmically. For 
simplicity, edge-on circular orbits with 10 random initial 
orbital phases were used. From analysing simulated light 
curves, the detectability of exomoons increases significantly 
with the moon's mass and decreases with increasing mass 
of planets. Moreover, the correlation coefficient of systems 
with the same planet mass or moon mass decreases with 
the planet semi-major- axis which can be explained by the 
correlation between the detectability and intensity of TDV 
signals which are weaker than TTV signals. Exomoon peri- 
ods in resonance with the planetary orbital period may pre- 
vent detection due to the constant observed planet orbital 
phase. The effects of intrinsic stellar variation (red noise) 
of M-dwarf reduce the detectability by 0.1. For simulation 
with red noise of system with planet masses less than around 
25 Me, 25-50% of simulated systems with 8-10 Me moon 
have correlations high enough to confirm the presence of an 
exomoon. 

In conclusion, the detectability of habitable exomoons is 
ultimately determined by the ability of Kepler to detect the 
TDV signal as this is typically weaker than the TTV signal. 
Resonance between planet period and moon period can pre- 
vent detection of some exomoon configurations. Planets with 
mass less than around 25 Me and star-planet orbital sepa- 
ration more than 0.3 AU should be the best candidates to 
detect exomoons in an M-dwarf system. Finally, exomoons 
in the habitable zone of an M-dwarf system can be detected 
by 3-years Kepler's data and the number of detected extra- 
solar moons should increase with the number of detected 
transiting exoplanets in the near future. 



i.o 
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• Iog(M ffl /M e ) = 1.0 

h- log(M m /iky = 0,9 
■ a log(M m /M m ) = 0.8 
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Figure 7. Median correlation as a function of planet of mass 
log(Mp/15M e ) equal to 1.0 M e (Blue circle with dot line), 0.9 
Me (Red square with dashed line) and 0.8 Me (Green triangle 
with dashed-dot line) exomoon orbiting around a M-dwarf with 
planet period and moon period are 89.35 and 2.24 days, respec- 
tively. Thick lines show the median correlations of the systems 
without stellar variability. The forward diagonal hatch region, 
shaded area and backward diagonal hatch region represent the 
25th to 75th percentile region of systems with 1.0, 0.9 and 0.8 
Me planets, respectively. 
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